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Abstract. In this paper we modified the Navier-Stokes equations by adding 
a higher order artificial viscosity term to the conventional system. We first 
showr that the solution of the regularized system converges strongly to the 
solution of the conventional system as the regularization parameter goes to 
zero, for each dimension d < 4. Then we show that the use of this artificial 
viscosity term leads to truncated the number of degrees of freedom in the 
long-time behavior of the solutions to these equations. This result suggests 
that the hyperviscous Navier-Stokes system is an interesting model for three- 
dimensional fluid turbulence. 



1. Introduction 

We regularize the Navier-Stokes equations by adding a higher-order viscosity 
term to the conventional system. In this paper we will restrict ourselves to periodic 
boundary conditions. 

^ -h e (-A)' We - vAue + (we.V) + Vp = f (x) , in X (0, oo) 

div — 0, in X (0, oo) , 
p{x + Lei,t) — p{x,t), u{x + Lei,t) = u{x,t) i — l,...,d te(0,oo) 
Ue (x, 0) = Ueo (x) , in ri, 

Where 51 (0, i)'' and (ei,...,ed) is the natural basis of R'^. Here e > is the 
artificial dissipation parameter and > is the kinematic viscosity of the fluid, 
I > 1. The function Ue is the velocity vector field, p is the pressure, and / is a given 
force field. For e = 0, the model is reduced to the Navier-Stokes system. 

In Lions [25], the existence and uniqueness of weak solutions of the modified 
Navier-Stokes equations were established for all Z > if / > ^j^, d is the space 
dimension. 

This type of regularization was proposed by Ladyzhenskaya [20] and Lions [26] 

who added the artificial hyperviscosity (—A) 2, I > 2 to the Navier-Stokes system. 

Mathematical model for such fluid motion play an important role in theoretical 
and computational studies of bipolar fluids [7] and in the regularized Navier-Stokes 
equations (see [T] [26] [28] and the references therein). Hyperviscosity has been 
widely used for numerical simulations of turbulence [U [3] [5] and in computer 
simulations for oceanic and atmospheric flows (see [4] [23]) or to control the Navier- 
Stokes equations [31j . 
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A well known example of such a result is the viscosity solution method for the 
Hamilton- Jacobi equations |26) . 

In this paper, we will study the effect of hyperviscosity on the Navier-Stokes 
turbulence. First, we show that the solutions of converge strongly to the 

corresponding solutions of the Navier-Stokes equations for d < 4. This result can 
extend to each domain D, with one finite size. 

In this result, we show that the conjecture of J. Lions [251 Remarque 8.2. SecII] 
is true, for d < 4. In addition, it is an extension of a result due to Lions 
(where only the weak convergence is proved). The results in this article can be 
seen as an improved version of the convergence results announced by Yuh-Roung 
and Sritharan [28l[29], in two different ways: On the one hand, we consider here a 
dimension d < 4, on the other hand the order viscosity term here is I > sup(|, ^j^). 

Next, we consider the system (|l.ip with / = 2 i.e. we modified the 3-D Navier- 
Stokes system by adding a fourth order artificial viscosity term (Laplacian square) 
and we show the existence of absorbing sets. This fact implies that the system 
{I = 2) possesses a global attractor Sl^. 

Finally, we obtain scale-invariant estimates on the Hausdorff and fractal dimen- 
sions of the global attractor 2le independent of e in terms of the Landau-Lifschitz 
theory [22 of the number of degrees of freedom in turbulent flow [HI [32] ■ In 
fact such an estimate that improves on the Landau-Lifschitz estimates has already 
been done by J. Avrin fTl in which hyperviscous terms are spectrally added to the 
Navier-Stokes equations. 

Thus we recover the improvement on the cubic power, i.e. get a bound propor- 
tional to G2 for p < 3. The latter should be a possibility, as the attractor results 
in [1] were not intended to be optimal in this direction. We would then represent 
an overlapping result that is new as far as we know, although readers familiar with 
the attractor techniques used may anticipate that such a result is possible in the 
hyperviscous case given the existing results in [I and the expected improvement in 
the Sobolev-space estimates in the fixed uniform hyperviscous case at hand. 

In Section 2, we present the relevant mathematical framework for the paper. In 
Section 3, we show the convergence of the system (jl.ll) to the conventional Navier- 
Stokes equations. In Section 4, we consider the hyperviscous system (l — 2), we 
show the existence of a global attractor. In Section 5, we estimate the dimension 
of the attractor. Finally, we provide in Section 6, explicit upper bounds for the 
dimension of the global attractor of the modified Navier-Stokes in terms of the 
relevant physical parameters. 

2. Notations and preliminaries 

In this section we introduce notations and the definitions of standard functional 
spaces that will be used throughout the paper. We denote by if™ (fi), the Sobolev 
space of L-periodic functions. These spaces are endowed with the inner product 

l/3|<™ 

and the norm 

ll"IL= E {\\D^u\\l 

l/3|<m ^ ' 

H-"' (fl) Denote the dual space of i?" (fl). 

We denote by 77™ (fi) the subspace of i/'" (il) with, zero average 
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i/™ = {u e iJ" (O) ; Ju{x)dx^O}. 

n 

For m = 0, we have ff™ (fi) = P (n). 

• We introduce the following solenoidal subspaces Vs, s E which are 
important to our analysis 

the spaces Vs as completions of smooth, divergence-free, periodic, zero-average 
functions with respect to the norms. 

Vq {fl) = {ueL^ (n) , divu = 0, u.n -u.n |s.+3 , « = 1, 2, 3}; 
Vi {n) = {ueH^ (fl) , divu = 0, 70U |e.= |e,+3 , « = 1, 2, 3}. 

V2 (fl) ^ {u e (il) ,divu = 0,7ou |s,= 70" |s,+3:7i" |s.= -7i" Isi+s-* = li2,3}, 

see [351 Chapter III, Section 2]. We refer the reader to R.Temam 33] for details 
on these spaces. Here the faces of fl are numbered as 

J:^^^nn {x., ^ O} and S.+a ^^nn{x^^ L}, i^l, 2, 3. 

Here 70, 71 are the trace operators and n is the unit outward normal on 

• The space Vq is endowed with the inner product ('u,u)^2(j2) and norm 

1 

\\U\\ = {U,u)l,^^y 

• Vi Is the Hilbert space with the norm — ll^lly ■ The norm induced by 

(17) and the norm ||Vu|| are equivalent in Vi. 

• V2 Is the Hilbert space with the norm ||u||2 = ll'^llv'j- Iii ^2 the norm 
induced by (fi) is equivalent to the norm ||Au||. 

Denote the dual space of Vs . We denote by A the Stokes operator 

Au = -Au for u e -D (A) . 

We recall that the operator A is a closed positive self-adjoint unbounded operator, 
with D (A) = {u e Vq, Au E Vq}. We have in fact, 

D{A) = in)r]Vo = V2. 

The eigenvalues of A are {Aj}^-~ij", < Ai < A2 < ...and the corresponding or- 

thonormal set of eigenfunctions {wj }^^^ is complete in Vq 

Awj = XjWj, Wj € D{Ai),\/j. 

The spectral theory of A allows us to define the powers A'- of A for Z > 1, 

A^ is an unbounded self-adjoint operator in Vq with a domain D{A^) dense in 

V2 <zVq. We set here 

A^u = (-A)' u for u e £) {A^) = V21 n Vq. 

The space D (A') is endowed with the scalar product and the norm 

("'f')D(A') = {A^U.A^v) , = {(U,w)z5(^i)}i 

Let us now define the trilinear form 6(., ., .) associated with the inertia terms 

b{u,v,w)— I Ui-—^Wjdx. 
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The continuity property of the trihnear form enables us to define (using Riesz 
representation Theorem) a bihnear continuous operator B (u,w); V2 XV2 ^ V2 will 
be defined by 

{B{u,v),w) ^b{u,v,w), Ww eV2. (2.f) 
Recall that for u satisfying V.m = we have 

b (u, u,u) — and b (u, ti, w) = —b (u, w, v) . (2-2) 

Hereafter, for i G N, will denote a dimensionless scale invariant positive constant 
which might depend on the shape of the domain. Similarly, the trilinear form 
b{u,v,w) satisfies the well-known inequalities (see, for instance, |30[ Lemma 61.1] 
and [81 (33!) 

\b{u,v,u)\ < ci llwllf llwll;^ for all u,vGV. (2.3) 

The trilinear form b{., ., .) is continuous on H"'^ (n) x (f^) x H"'-- (n), 

rrii > 

3 

|6(m,u,w)| < C2 ||u||„^J|?;||„^^_^i ||w||,„^ , + ni2 + mi > - (2.4) 

see [21) . We recall some well known inequalities that we will be using in what 
follows. 

Agmon inequality (see, e.g., [5]) 

||u|L < C3 hllf \\Au\\^ for all u £ V2. (2.5) 

Young's inequality 

ab<^aP + ^b\a,b,a>0, p > 1, g = (2.6) 

Poincare inequality 

Ai < ||A5uf for all uGV. (2.7) 

To prove uniform bounds on different norms we use the uniform Gronwall Lemma 
for proof see [321 Lemma III 1.1]. 

Lemma 2.1. (The Uniform Gronwall Lemma) Let g, h, y be three positive locally 
integrable functions on (to, + 00) which satisfy 

^ gy + h for t > tQ and 

t+r rt-{-r ft+r 

g{s)ds<ai, I h{s)ds<a2, / y {s) ds < for t > to, 
where ai, 02, 03 are positive constants. Then 

y{t + r) < + 02^ exp (ai) for t > to. 

Denoting by G the dimensionless Grashoff number 10 , G = in 3D, (see 

e.g. [TJ [TTJ [33]). Tthis number measures the relative strength of the forcing and 
viscosity. 
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3. Strong convergence for the hyperviscous system 

In this Section, we give a new Theorem which insures the strong convergence of 
the solutions of the system (|l.ip to the corresponding solutions of the Navier-Stokes 
equations for d < 4. This result can extend to each domain with one finite size. 
Moreover, we show that £ C (0, T; Vq). 

Using the operators defined above, we can write the modified system ()1.1|1 in the 
evolution form 

— + eA^u, + uAu^ + B (u„ u.) ^ f (x) , in il x (0, oo) (3.1) 
dt 

UeO {x) = u^o , in Q,. (3.2) 

The existence and uniqueness results for initial value problem (jl.ip can be found 
in [25], m Chap.l, Remarque 6.11]. 

The following theorem collects the main result in this work 

Theorem 3.1. For I > d is the space dimension, for e > fixed, f £ L"^ (0, T; V^) 
and € Vq be given. There exists a unique weak solution of (| j. which satisfies 

Ue e (0, T; Vi) n (0, T; Vb) , VT > 0. 

Notice that the conventional Navier-Stokes system can be written in the evolu- 
tion form 

du 

— + uAu + B{u,u) = f (x) , in ^1 X (0, oo) (3.3) 
u (0) = uq , in n. (3.4) 

Theorem 3.2. For d < i, for f e L^(0; T; Vq) and uq S Vq be given. There exists 
a weak solution of ((23)-([^ which satisfies u G L°°(0;T;Fo) n L'^[Q;T;Vi), for 
T > 0. For d — 2, u is unique (J. Lions |25| j. 

We will establish various estimates uniform in s for the solutions of the modified 
Navier Stokes. These bounds will be used to establish the limit of these solutions 
to the conventional Navier Stokes equations. 

Proposition 3.3. For c? < 4 and for e > fixed, f & L'^ (0, T; Vq) and u^q € Vq. 

The weak solution (t) of the modified Navier-Stokes equations satisfy 

i) Ue is uniformly bounded in L°° (0,T; Vq), 

ii) Ue is uniformly bounded in {0,T;Vi). 

We need the following Lemma proved in R. Temam j33l Lemma 4.1.ChIII,Sec4]. 
Lemma 3.4. The form b is trilinear continuous on V x V x Vs if s > ^ and 

\\b{u,v,w)\\ < C4||u||||t;||i||w||5. 
Applying Lemma 13.41 we obtain 

Lemma 3.5. Let [t) be a weak solution of the modified Navier-Stokes system. 
Then B (u^) belongs to (Q, T; VJ') forl>^. 
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Proof. By the definition of the operator B and the above Lemma, we get 

\{B{u{t),v))\ = \b{u{t),u{t),v)\ 

<c4u{t)\\\\u{t)Uv\\v>, WveVi. 

We set 

B{u{t)) = B{u{t) ,u{t)), 

thus 

\\B{u{t))\\v,, <c4u{t)\\\\u{t)\\i ioTO<t<T. 

□ 

Lemma 3.6. If f £ (0,T; ¥{), then, for any solution u^ (t) of problem the 
du 

time derivative is uniformly hounded in (OjT; V^'). 

Proof. Due to Lemma [?31 i3 (u^) belongs to (0, T; V/), since f—eA^u^—vAu^ belongs 

du 

to £2 (0, T; VI), this implies that belongs to (0, T; □ 

Lemma 3.7. is almost everywhere equal to a continuous function from [0,T] to 
the space Vq. 

du 

Proof. Since u^ G (0,T;yi) n L°° (0,T;Vb) and -j^ ^ (fi,T;V{), the weak 

continuity in Vb is a direct consequence of [321 Lemma 1.4.ChIII,Secl]. 

Similarly, it follows that u^ (0) converges to u (0) in Vq , and since u^o converges 
to uq in V/ , we conclude that u{0) = uq- D 

Now we prove the strong convergence. It follows from ii) of Proposition 13.31 and 
from Lemma |3.6[ that easily 

du^ 
~dt 



u,^ G X G L2(0,T; Vi), G (0,r; V^')} with bounds independent of 



En- Hence 



du 



(i) u,^ ^ u in L^{0,T;Vi) weakly; and {ii) ^ "T -^^(0. weakly; 

These two properties allow us to establish the strong convergence result. 

The proof of the following theorem can be found in R. Temam [331 Theorem 2.1, 
Chapter III, Sec 2]. 

du 

Theorem 3.8. The injection of X ^ {u e L'^ {0 , T ; Vi) , e (0,r;V,')} into 
y — {u ^ L^{0, T; Vo)} is compact. 

By virtue of the above estimates and the compactness Theorem l3.8l We can now 
state our first result. 

Theorem 3.9. For I > sup(|,^i^) and for d < 4, the weak solution u^ of the 
modified Navier-Stokes equations <\1.1]\ given by Theorem \3.1\ converges strongly in 
1^2(0, T; Vo) as e ^ to u a weak solution of the system l^3.3^-(3^. 



Proof. Theorem 13. II and Lemma [3.41 are satisfied for / > sup(^, ^-^). We use part 
ii) of Proposition 13.31 and Lemma 13.61 we can deduce that the weak solutions 

du 

e X ={u,^ G L^{0,T;Vi), G (0, T; V,')}. Hence, the compactness 

Theorem 13.81 implies the strong convergence in £^(0, T; Vq). □ 
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The following proposition is a consequence of Proposition 13.31 
ion 3.10. Vw e L^{0,T;Vi) 

T.dUe^{t) „T,du{t) 



Proposition 3.10. Vw G (o, T; Vi), £ (q^ y. 



h) lim (Vu£„ (t) , Vw (t)) dt = {Vu {t) , Vw (t)) dt, 
c) lim fl' b (Me„ (t) , Ue„ (i) , u) (t)) dt ^ 6 (u (t) ,u(t),w (t)) dt. 



Let us now establish the limit of the equations (I3.ip as e„ — 0. Taking the inner 
product of p.ip with a test function G 25(0, T; 2?(A5)) then integrate by parts 
and using the convergence Proposition 13. lOl we can pass to the limit as e„ — > 0, we 

T T T T 

get - /g (m, Lp') dt + vj^ (Vw, Vip) dt + J^ b {u, u, ip) dt = (/, i^) dt. 
Here the term e„ J^lA^Ui,^ (t) , A^ip [t))dt goes to as e„ — 0. 
Since the weak solution u^^^ £ (0, T; Vi) with bound uniform in £„ and we get 



En /(f {A'^Ue^.A^ip) dt<en \{ue„, A'' ip) \ dt < CEn- 

Since u & (0, T; Vl) n L°° (0, T; Vb), we can conclude that u is indeed a weak 
solution for the conventional Navier- Stokes equations. 

4. The hyperviscous Navier-Stokes system and attractors 

Now, we consider modifications of the 3D Navier-Stokes system by adding a 
fourth order artificial viscosity term (Laplacian square) depending on a small pa- 
rameter £ to the conventional system. 



i-T 



— h eA^u^ + vAue + B (u^, u^) ~ f {x) , in x (0, oo) 

dt (A\\ 
div Ue — Q, in fi X (0, oo) , Ue {x, 0) — Ueo (x) , in il, ^ ' ' 

p{x + Lei, t) = p{x, t), u{x -\- Lei, t) — u{x, t) i ^ 1,2,3. t G (0, oo) 

where ft = (0, L)^ . In this section we will show the existence of the compact global 
attractor St^ associated with the semigroup (<) generated by the problem (|4.1I) . 
(For the theory of global attractors see 0, [5], [E], [IB], ED, ED], [12] ■)• 

For e = weak solutions of problem are known to exist by a basic result by J. 
Leray from 1934 [24], only the uniqueness of weak solutions remains as an open 
problem. Then the known theory of global attractors of infinite dimensional dy- 
namical systems is not applicable to the 3D Navier-Stokes system. 

The theory of trajectory attractors for evolution partial differential equations was 
developed in [30] , which the uniqueness theorem of solutions of the corresponding 
initial- value problem is not proved yet, e.g. for the 3D Navier-Stokes system (see, 
for instance, [14i I30j ). Such trajectory attractor is a classical global attractor but 
in the space of weak solutions. 

The problem of upper semicontinuity of global attractors for the 2D with periodic 
boundary conditions was discussed by Yuh-Roung Ou and S. S. Sritharan in [28) . 
For related results which use the theory has been introduced by Foias, Sell, and 
Temam in [12l |32] to show that the system (|l.ip possesses an inertial manifold (see 
[D [291 [32]). 

The existence and uniqueness results for initial value problem (|4.1|) are conse- 
quence of Theorem [Hm for I = 2 and d — 3. 
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Theorem 4.1. Let C M^, and let f e L'^ (0,r; Vj) and Ueo e Vq be given. Then 
there exists a unique weak solution of ( |^. J[ ) which satisfies 

£ C ([0, T] ; Vo) nL^ (0, T; V2) , VT > 0. T/ien as e — > 0, the solution converges 
to a weak solution of the Navier-Stokes equations. 

Now, we show that the semigroup Ss (t) has an absorbing baU in Vo and an 
absorbing baU in Vi. Then we show that (t) admits a compact attractor in Vq 
for each £ > 0. 

We take the inner product of (|4.ip with u^, we obtain the energy equahty 
^ \\u,f + 2s \\Au,f + 2v ||Vu,f = 2 (/, w,) . 

Here we have used the fact that h (u^, u^, u^) — 0. By applying Young's inequality 
and the Poincare Lemma, we get 

4 \Wf + 2£ \\Au£ + V ||Vu,f < (4.2) 

we drop the term 2e we obtain 

by integrating the above inequality from to t,we get 

\W (i)f < ll"eof e-''^^* + Po (1 - e-''^^*) , t > 0, (4.3) 

where po = —7— ||/||- Hence for any ball Br^ — {u^o G Vo; ||Meo|| < ^0} there is a 
ball i? (0, bo) in Vb centered at origin with radius 60 > Po (^^o > (^o) such that 

^logt^ 

The ball is said to be absorbing and invariant under the action of S^(t). 
Taking the limit in (14.31) we get, 

lim sup ||we(t)|| < po. (4.5) 

t—^OO 

We integrate (|4.2D from ^ to i + r, we obtain for w^q ^ 

r lke||?ds< i(^^^ + h,(i)f), Vr>0, Vi>io(i3flo)- (4-6) 



S, {t)BR„ C Br„ for t > to {Br„ ) = — log ,° 2 • (4-4) 



With the use of (|4.5I) we conclude that 



from which we obtain 



lim sup / \\u4lds < ll/f + M!, (4.7) 



1 



lim sup - / ||w,||2ds< (4.8) 

this verifies that the left-hand side is finite. 

To show that the semigroup Seit) has an absorbing set in Vi^ we consider the 
strong solutions and take the inner product of (|4.1[) with Au^, we obtain 

\j^\\A^u,f+e\\A^u,f + y\\Au,f = -b{u,,u,, An,) + {f,Au,). (4.9) 



EFFECT OF HYPERVISCOSITY ON TURBULENCE. 



9 



By applying Young's inequality, we get 

if,Au,)<\\f\\\\Au,\ 



4 I' 



By using the Agmon's inequality (|2.5p and Young's inequality we can estimate the 
last term in the left-hand side of (14.91) as follows 



\b{u,,u,,Au,)\ < \K\\^\K\\^\\Au,\\ 
<C4||u,||f \\Au,\\i 
<!i\\Au,f+c,\K\\l 

Hence we obtain from (|4.9p 

^ + 2e\\Aiu,r + y \\Au£ < - ll/f + 2c5 . 
at V 

Dropping the positive terms associated with e we have 

+ V \\A,u£ < "Wl. + (4.10) 
at V 

we apply the uniform Gronwall Lemma to (I4.10p with 

5 = 2c4||ti,||l, h^'^-^,y=\Wt. 

For n = 3, m = 2 and = i, in [26| Formula (6.167)], we get = 6 wich means 
e (0, T; Vi) then e L"* (0, T; Vi) , thus 

Thanks to (|4.3p - (|4.7p we estimate the quantities ai, 02, a3 in Gronwall Lemma by 

2r\\ff r ll/f ll/f 
ai = 20404, 02 = , 03 h ^r-TT- 

Then we obtain 

\\ue {t)\\l < (y +a2)exp(ai) = Rj for i > io, as in (lO) . 

Hence, for any ball , there exists a ball -65^ , in Vi centered at origin with radius 
i?i > (5i > pi such that 

1 E? — 

Se{t)BR, C for t > ti (BflJ = to (5i?o) + 1 + ^ log ^. 

I^Ai Oi — Pi 

The ball Bg-^ is said to be absorbing and invariant for the semigroup S'e(t). 

Furthermore, if B is any bounded set of Vb, then Se{t)B C for t>ti {B, Rq), 
this shows the existence of an absorbing set in Vi. Since the embedding of Vi in 
Vq is compact, we deduce that Se{t) maps a bounded set in Vq into a compact set 
in Vq. In addition, the operators S^(t) are uniformly compact for t > ti {B,Ro). 
That is, 

is relatively compact in Vq. 
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Due to a the standard procedure (cf., for example, |32l Theorem 1. 1.1] for de- 
tails) , one can prove that there is a global attractor a compact attractor 2te for the 
operators S'e(t) for £ > 0, 

Note that the global attractor Slj must be contained in the absorbing balls Vq 
and Vi 

2te = nt,>oUt>uBs^ it) C Bs, n Bs, . (4.11) 
Notice that all the above bounds are independent of e. 



5. Estimates of Dimensions of the Global Attractor 

Our aim in this section is to study the finite dimensionality of the global attrac- 
tor. In the first part we will prove the differentiability property of (t) and in the 
second part we will provide estimates of the fractal and Hausdorff dimensions of 
their global attractors Sl^. 

Using the trace formula [32l Chapters V and VI] , we estimate the Hausdorff and 
the fractal dimensions of the global attractor 2le in V. 

For a solution (t) = (t) u^o, t > 0, lying on the attractor u^o G 21^, we see 
from (|4.ip that the linearized fiow around is given by the equation 

U^+eA^Ue + iyAU, + B{ue,Ue)+B{Ue,u,)^0, inV , . 

U, (0) = in V. ^ ' 

We show the differentiability of the semigroup Se with respect to the initial data 
in the space V. 

Theorem 5.1. For any t > 0, the function u^q — > (t) — (t) Weo 'is Frechet 
differentiable on the attractor 21^ . Its differential is the linear operator 

D {S, (t) u,o) - L (t, u,o) -.CeV ^U,{t) eV,te[0,T], 

where Ug (t) is the solution of . 
Proof. Let 

7^{t) = Ve (t) - Ue (t) - Ue {t) , C/e (0) = ^ = "eO " UsO- 

Clearly, r] satisfies 

77t + eA^-q + vArj + B{rj, v,) + B{ve,ri) - B{w,,w,) = 0, ry(0) = 

where = —Ug. Taking the inner product of the last equation with 77 and using 
the identity Biv^, r],r]) = we obtain 




+ 2e \\Ar)f + 2v \\vif^ = 2&(7?, w„ 77) - 26K, u;„ rf). (5.2) 



By (|2.3p the first term in the right-hand side of (|5.2p has the estimate 
|26(7?,z;„77)|<2ci||,7p||r;||f 

<2cii?ih|i^h||f 

<-^hll + — hill- 
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Employing the inequalities (|2.3|) we estimate the second term in the right hand side 
of (15.211 as follows 



2b{we,We,r]) < 2ci ||?7||^ \\w, 



2c? ,, ,,4 ly 

<^IKIlt+2""-- 



elll 

2 



Hence, we obtain from (E 



d\\vf , o II 4 i|2 , 3;^ ||2 ^ cjRj II ||2 , 2c? 4 
— ^ + 2e||Ar7|| + — II77II1 < h|| + — 

we drop the positive terms 2e HAryU^ and ^ ||77||J we get 

^<^||.f + ^IKIlt. (5.3) 
From the classical Gronwall Lemma (see |33])j (|5.3p gives 

„2 /-t ^4 d4 



,2 2cf 

r < — i 



/ l|u'e|ltexp( / ^-5J-rfr)rfs 



V 

Thus 



<cJ IKII^ds, Co = ^exp(^i^ 
Jo V V 



\2h{we,v,,We)\<2cx\\v,\\^\\wM\\we\\'- 
c\R\ ,, ,,2 , 3i. 



<^Kir + -iKii^. 



(5.4) 



The difference 

We (t) = Ve (t) - (t) = S'e (i) V^o ~ Se{t)UeO 

satisfies the equation 

^ + eA^We + vAwe + B{we,ve) + B{ve,We) B{we,We) =0 (5.5) 
at 

and 

Taking the inner product of the last equation with We,we obtain 

^ \\w,f + 2e \\Aw,f + 2u \\w,\\l = 2b{w„w„ v,). (5.6) 
By using inequalities (j2.3p and Young's inequality we obtain 



,3 ' 4-""-" 

Substituting the above result into (|5.6p . we obtain 

I + 2e WAweW' + ^ ||«;,||? < 1^ ||^,|^ (5.7) 

We drop the positive terms 2e \\Awe\f' and ^ \\we\\i to obtain the following differ- 
ential inequality 

||Kf <^Kf . (5.8) 
Using the classical Gronwall Lemma we deduce from (|5.8I) that 

lkeir<lke(0)irexp(^^). (5.9) 
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From (I5.9P we deduce that 

\K{t)~v,it)\\ldt<C,\\u,o~v,of;C, = ^rexp(^^), (5.10) 
with (|5.4p we conclude that 

hf <CoC2||7.,o-Weo||', 

then we deduce from (|5.4|) and (|5.10p that 

Ik/f < C2 Ih, (0)|i\ where C2 - CqC? (5.11) 

this shows that 

The differentiabihty of (t) is proved. □ 



From Theorem 15. li the function (f) is Frechet differentiable on 2le for t > 0. 
For ^ G Vb, there exists a unique solution of (|5.ip satisfies 

{/eeC([o,r];Fo)nL2(o,r;T/2) vr>o. 



With the differentiability ensured in Theorem 14.11 we can then define a linear 
map L {t; u^q) : ^ G Vq — > J/^ (t) £ Vq where Ug is the solution of (15. ip . 

We can apply the trace formula (see [5] and \W2\ Section V. 3]) to find a bound 
on the dimension of the global attractor 21^. We consider the trace TrF' (u^) of 
the linear operator F' (u^) and for to G N, the number 

1 /■* 

= lim sup sup sup - / TrF' {S^ [t) u^o) o Q^{t) dr 
l?il<i 

where Qm (t) = Qm (t, Weo; ^ii ?m) is the orthogonal projector in Vo onto the 
space spanned by (r) , C/™ (r). where {t) = L{t,Uso) .^j, j — l,...,m, 
t > 0, are to solutions of ()5.ip . corresponding to ^ = Ci,...,Cm G ^i- Let 
ipj{T), j — 1,...,TO, r > 0, be an orthonormal basis of for Qm (r) Vq =span 
{[/i(r),...,C/r(r)}, ^,(t) € 1/1 for J = 1,...,TO, since {t) , ...,U^ (r) e Fi, 
T G R+. 

From the general result in [32l Section V.3.41], we have that if qm < for 
some m € N then the global attractor has finite Hausdorff and fractal dimensions 
estimated respectively as 

dimn (2te) < to, (5.12) 
dinii. (2le) < to(1 + max -[p^). (5.13) 

l<i<m-l \\q,n\\ 

Then we have 

TrF' {S, (r) u,o) o Qm (r) = E.'li (^ri^' ("s (r)) o (r) (r) , cp, (r)) 

m 

= E (^^'(ti, (r))^,(r),^,(r)). 

i=i 

Recall that (., .) denoting the scalar product in Vq, we write using (|2.ip and (|2.2p 
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Tr{F' {u, {t))lPj (r) , tp^ (r)) = {s^'^'Pj - ^^Vj - B{ipj,u,) - B{u,,ifj),ipj) 

2 1 

= EC-ell^V'jII -^WA'^V'jf -b{ue,(Pj,(Pj) -b{ipj,'. 

thus 

Tr {F' (u, (r)) (r) , (r)) = E;'Li(-£ 11^. II2 " ^ Hv. II? " ^ (^., ^j))- (5.14) 
We estimate the nonlinear term as follows 

m m P 3 Qy^ 

I Y.b{Vj,u,ipj) HJ2 E ^jzjr^ (x) (fijkdx \ (5.15) 

j"=l j=l Jo i,k=l CXi 

whence for almost every x G Q we have 

m 3 Qy^ 

I E E "^JiTd^ (^) <^jfc^2; |< ||w||i IIpII 

J — 1 ^Xi 



where 



and 



k(^)iii = ( E muAx)f)^ 



m 3 

EE 

i=ii=i 



p{x)^Y.Y.{^o^{^)f ■ (5-16) 



Therefore 

m r 

\Y.b{Vj,Ue,^Pj)\< / p{x)\\ue{x)\\-^^dx. (5.17) 
J=l Jo 

Now we recall the generalized form of the Lieb-Thirring inequality in dimension 

three and m = / as developed in [32 , Theorem A4.1] 

Theorem 5.2. (The Lieb-Thirring inequality). Let (pj, 1 < j < N be a finite 
family ofVi wich is orthonormal in L^{fl) and set, for every x G fl, 

p(^) = E \\ivA^))f 

3 = 1 

Then there exists a constant k, independent of the family ipj and of N such that 

( / p{xf^-'dxf(''-^y^ <nY. I a{^„^,). (5.18) 
Jn j=i Jn 

for all q G max{{l, S/21), (1 + 3/21)} and where k depends on I, p, and q, and on 
the shape (but not the size) offl. 

The quadratic form we will use is 

a{v, u) = {A^i, u) = (A'/2w, A^/^u) (5.19) 

so that the order of our quadratic form is I. 

Kolmogorov's mean rate of dissipation of energy in turbulent flow (see e.g. |lll 
VI. (3.20)]) is defined as 

3 1 /■* 2 

e = Af* lim sup sup - / ||ue(r)||^dT (5.20) 
t->oo ueoeae ^ Jo 



14 



ABDELHAFID YOUNSI. 



the maximal mean rate of dissipation of energy on the attractor, which is finite 
thanks to (j^S]) . 

Using (|4.8I) we can estimate the energy dissipation flux e by 




In order to make the dimension estimate more exphcit, we can estimate the energy 
dissipation flux e in terms of G by 

e < Xlv'^G^. (5.22) 



6. Numbers of degrees of freedom in turbulent flows 

In this Section, we estimate the effects of hyperviscosity on the turbulent flow. 
An argument from the classical theory of turbulence (see, L. Landau and Lifshitz 
|22| ) suggests that there arc flnitcly many degrees of freedom in turbulent flows. 
Heuristic physical arguments are used to justify this assertion and to provide an 
estimate for this number of degrees of freedom by dividing a typical length scale 

of the flow, /q = ^ , by the Kolmogorov dissipation length scale 1^ i.e. = ^ 
where e is Kolmogorov's mean rate of dissipation of energy in turbulent flow and 
taking the third power in 3D. 

We will express our primary attractor results in terms of the Kolmogorov length- 
scale and the Landau-Lifschitz estimates [52] of the number of degrees of freedom 
in turbulent flow [TTl [32j and we can easily observe such compatibility that exists 
between these estimates and the number of degrees of freedom in turbulence (see 
also [22]). Such estimates will give us useful information about the capability of 
(|4.1|) to approximate Navier-Stokes equations dynamics. We will show that the 
corresponding number of degrees of freedom is proportional to the dimension of the 
global attractor. 

By Holder's inequality the right hand side of (|5.17l) can be estimated as follow 
\\u,{x)\\^p{x)dx<\\p{x)\\^.^^^^ \\A^u,{x) ll^j^^^ (6.1) 
Applying Young's inequality with 



we obtain 

/ \\ue {x)\\^p{x)dx < ^ ||p(a;)||% +C5 || A^u^ (x) || S , cs = ^i^) 
Jq 2k L3(o) L4(f2) 7 6k 

Using (|6.3|) we can majorize TrF' {u^ (r)) o (r) as follows 



4,7e._3 

4 

(6.3) 



TrF' (u, (r)) o Q„ (r) < -^E ||^, {x)\\t - 11^. Mll^ + 11/^(^)11^.0. 

(6.4) 

+c,\\A-^u,ix)\\\^^^. (6.5) 
Applying the Lieb-Thirring inequality (|5.18l ) we obtain 
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TrF'{u,{T))oQ,^[r) < --^ \W,{x)\\l--Y. 11^, (r) ||^ + C5 || A-^u,{x) ||% . 

Z j—i L j—i l^i (it) 

The Sobolev embedding V2 C Vi the Sobolev inequahties on 57 in terms of 

II (x) ||l< Cg II ifij (x) II2 

we get 

TrF' {u, (r)) o Q„, (r) < + — ) E ll^^, (^)ll + C5 II A'^u, [x) f . . (6.6) 
then, 

TrF' {ue (t))oQ„ (t) < -C7I: llv^j {x)\\\+c^ \\ A'-u, {x) \\\^^^ , where cj = - + — . 

(6.7) ' 

Note that in the 3D case we have > csL~'^j's for some positive universal 
constant (see, for example [351 Lemma VI 2.1]). Therefore, 

E \\ipj {x)\\l > Ai + ... + A,„ > cgAima . (6.8) 

For the term || Azu^ (x) || * 7 , we have by Holder's inequality that 

L 4 (o) 

II (x) ||J|(^^< cio II A^Me (x) II' with CIO = (6.9) 



Taking into account (|6.4p then yields 

Tri^' (u, (t)) o Q„ (r) dr < -cyCgAimi + C5C10 || A^u, (x) \\i . (6.10) 

By Holder's inequality we have 

1 /■* 1 7 1 /■* 1 1 

hm sup sup - / II A'^u^ {''':x) |1 * dr < hm sup ( sup - / || A'^u^ {T:X) |P dr)'' 

t— fOOM^oGate t Jo t— »CX3 u^o esie i Jo 

(6.11) 

On the other hand, using ()5.20p we have 

1 /■* 1 7 £7 

lim sup sup - / \\ A^Ui; {t,x) \\i dr < . (6-12) 

For u^Q G 2te, we can estimate the quantities qm (t), Qm 



where 



Qm ~ lim sup Qm (t) < — Kim^ + K2, (6.13) 



e 7 

Ki = C7C9A1 and K2 = C5Cio(^— )s. (6.14) 



Therefore, if m' £ N is defined by 



/ 1 ^/-^'^2n3 . 2C5C1O .3 21 /^ir-N 

m - 1 < ( )5 = ( 37 )5e4o < to', (6.15) 

'^1 CyCgAj^l^i 
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3 1 _i 

Setting — ( — )3 the dissipation length scale, and = -^i ^ the macroscopical 
length by setting. Then we can rewrite (|6.15p in the form 

m' - 1 < c,,(^)TO < m' (6.16) 

where 

c. = (^^)«(.-)Af. (6.17) 
Thus, we have proved the following Proposition 

Proposition 6.1. The Hausdorff and fractal dimensions of the global attractor 21^ 
of the regularized 3D Navier-Stokes l \4-lh dimp (21^) and dimn (Sl^) respectively, 
satisfy 

dimH(2t,) < dini;^(2le) < c,A^. (6.18) 

The exponent on js. is significantly less than the Landau-Lifschitz predicted value 
of 3, less than the results in ^ for the 3D Camassa-Holm equations, or simply NS-a 
model and less than the Avrin exponent (for a = I = 2) [1, Theorem 1]. 

This, in a sense, suggests that in the absence of boundary effects (e.g., in the 
case of periodic boundary conditions) the modified 3D Navier-Stokes represent, 
very well, the averaged equation of motion of turbulent flows. 

11/11 

Since the Grashoff number G — 3- in 3D, (see e.g. pQ[TTl[33]) is an upper 

bound for (|^)^, expressing the above estimates in terms of G is straightforward. 
The above Proposition becomes 

Proposition 6.2. The Hausdorff and fractal dimensions of the global attractor Sl^ 
of the regularized 3D Navier-Stokes ( |/^._?p , dim^? (2le) and dim^f (2le) respectively, 
satisfy 

diniH (2le) < dimj. (2le) < c,,G5n. (6.19) 

Thus we recover the improvement on the cubic power, i.e. get a bound propor- 
tional to G2 for p < 3, in (I6.19P P = fg. This improvement suggesting to very good 
agreement with the conventional theory of turbulence. 

For a = I ^ 2, motivated by the Chapman-Enskog expansion, we recover (|6.19p . 
This result can be seen as an improved version of the results announced by Joel 
Avrin [1] Theorem 2]. 

This upper bound is much smaller than what one would expect for three-dimensional 
models, i.e. ^tf. This improves significantly on previous bounds have demon- 
strated that hyperviscosity can have profound effects on the number of degree 
freedom. The modifying effects are well understood, which makes the use of hyper- 
viscosity an efficient tool for numerical studies and suggests that the regularized 
3D Navier-Stokes has a great potential to become a good sub-gridscale large-eddy 
simulation model of turbulence. The results obtained agree very well with those 
provided in numerical studies of turbulence(see, Refs.,[T],[5], [T3], [IS], PT]V 

The present results explain some fundamental differences between the theory use 
instead a hyper-viscous term to approximate Navier-Stokes equations and which 
hyperviscous terms are added spectrally to the standard incompressible Navier- 
Stokes equations [1]. It would be interesting to obtain estimates for ()l.ip in this 
context in 3D and to see how the estimates depend on ^ for / > |. 
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